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We study the scaling of decoherence of decoupled electron spin qubits due to hyperfine interaction. 
For a superposed state consisting of product states from a single Zeeman manifold, both Tf{n) and 
T 2 (n) are scale-free with respect to n and the number of basis states, m. For a superposed state 
made up of states from different Zeeman manifolds, both ? 2 (w) and Tf{n) are roughly inversely 
proportional to ydi. Our results can be extended to other decoherence mechanisms, including in the 
presence of dynamical decoupling, which allow meaningful discussions on the scalability of spin-based 
coherent solid state quantum technology. 


Introduction. —Large-scale quantum information pro¬ 
cessing (QIP) requires the generation, manipulation, 
and measurement of fully coherent superposed quantum 
states involving many qubits [1] . One of the key issues for 
QIP is how well such a many-qubit system can maintain 
its quantum coherence. This is also an important issue 
from the perspective of fundamental physics: it remains 
an intriguing question how a large number of microscopic 
quantum mechanical systems together behave classically 
as a macroscopic object. Again, decoherence is central 
to such quantum-to-classical transitions [2] . 

A confined single electron spin in a semiconductor 
quantum dot (QD) or a shallow donor is highly quan¬ 
tum coherent, and is an ideal candidate as a qubit [3-7]. 
At low temperatures, an isolated electron spin has an ex¬ 
ceedingly long longitudinal relaxation time [8-11], and a 
very long pure dephasing time after removing inhomoge¬ 
neous broadening [12-15]. It is now well understood that 
the main single-spin decoherence channel is through hy¬ 
perfine coupling to the environmental nuclear spins [13- 
17], and the effects of hyperfine interaction have also been 
investigated for coupled two-, three- and even more spin 
systems [18-26]. On the other hand, decoherence of a 
many-spin system is still unexplored. 

In this Letter we study hyperfine-induced decoherence 
of n {:s> 1) decoupled QD-confined electron spin qubits. 
Our goals are to clarify how fast a many-qubit superposed 
state loses its coherence, and how this collective decoher¬ 
ence scales with the number of qubits involved. In our 
study, a uniform magnetic field is applied, so that the 
Zeeman splitting is much larger than the nuclear-spin- 
induced inhomogeneous broadening (see Fig. 1). Conse¬ 
quently, the dominant single-spin decoherence channel is 
pure dephasing due to the nuclear spins [16, 17]. We ex¬ 
plore how this dephasing mechanism affects a many-spin- 
qubit state by examining a large number of superposed 
states in various forms. Our results from this broad- 
ranged exploration indicate a sublinear scaling behavior 
for dephasing rates in the short time limit, making the 
scale-up of a spin-based quantum computer a difficult but 
not intractable endeavor. 

Electron-nuclear spin hyperfine interaction. —We con¬ 
sider n decoupled electron spins in a finite uniform 
magnetic field, each confined (in a quantum dot, nom- 
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FIG. 1. The sketch about the energy distribution of manifolds 
{0 < k < n) for n electrons separately confined in n decoupled 
QDs. Each electron spin is coupled to local nuclear spins by 
hyperfine interaction that gives rise to local field as large as 
AS, so that the energy level for each manifold is broadened 
to a band with width ^/nAB. 


inally) and interacting with a local and uncorrelated 
nuclear-spin bath through hyperfine interaction. The to¬ 
tal Hamiltonian for this n-qubit system and the nuclear 
spin reservoirs is 

Htot = Hze + Hz„ + 

\ «=i 

where Ul is the electron Zeeman splitting, ujja is the nu¬ 
clear Zeeman splitting of the a-th nuclear spin in the 
j-th QD (from here on j will always be used to label the 
QDs), and Aja is the corresponding hyperfine coupling 
strength. The number of nuclear spins coupled to each 
electron spin, Nj, is generally large, in the order of 10® 
to 10® in GaAs QDs, and ~ 10® in natural Si QDs. 

With the electron spins isolated from each other, the 
total Hamiltonian is a sum of n fully independent single¬ 
spin decoherence Hamiltonians. The evolution operator 
for the n-qubit can thus be factored into a simple product 
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of operators for each individual qubits (before and after 
tracing over the local nuclear reservoirs). We present 
a brief recap of single-spin decoherence [16, 17] proper¬ 
ties in appendix A, and focus here on how we approach 
the multi-spin-qubit decoherence problem based on the 
results of the single-qubit case. We also note that in¬ 
homogeneous broadening and the narrowed-state free in¬ 
duction decay are statistically independent because of 
independence between longitudinal and transverse Over- 
hauser fields, as presented in appendix C. These two 
pure dephasing channels follow the same scaling law, i.e., 
T 2 *(n)/T 2 *(l) = r 2 (n)/T 2 (l). Thus in the following we 
will focus on the scaling analysis of T 2 (n). 

Multi-spin decoherence. —For an n-spin system in a 
finite uniform external magnetic field, the full Hilbert 
space is divided into n -|- 1 Zeeman subspaces, labeled by 
the expectation value of = fc — n/2, fc = 0, 1, 2, • • • , n. 
Each subspace consists of = n\/[k\{n — fc)!] degener¬ 
ate states (in the absence of nuclear field), which has k 
spins in the jl) state and n — k spins in the jl) state. 
The local random Overhauser field breaks this degener¬ 
acy and leads to a broadening of ^ y/nAB, as illustrated 
in Fig. 1. For our decoherence calculations, we use the 
spin product states \xr) = \ln^n-i ’ ”^i) the bases. 
Here Vj refers to the electron spin orientation along the 
^-direction in the jth QD for state ja;,.), and takes the 
value of 1 or 1 = —1 for notational simplicity . 

For a superposed state ja;) that contains more than 
one product state, decoherence emerges due to the non¬ 
stationary random phase differences among m product 
states |a;r.)’s: \x{t)) = exTp{—iBf.rir ^...i^t}\xr) 

with 1 (from now on we use m to repre¬ 

sent the number of product states contained in ja;), and 
the notation for the Overhauser field is defined in the 
appendix B). As a collective decoherence measure of ja;) 
caused by the inhomogeneous broadening [i.e., in the fol¬ 
lowing calculations we use only the longitudinal Over¬ 
hauser field B instead of the total Overhauser field B], 
we use fidelity defined as B{t) = Y^M[(a;|a;(t))(a;(t)|a;)] 
[see Eqs. (C2) and (C3)]. For ja:), we find 


\dr\'^ + 2^\ddkVM[cos{e,kt)] , (2) 

r=l i<ck 

where the phase difference is Oik = ^ ^ — 

Bfkjk ,k- Specifically, M[cos(6ikt)] is solely deter- 
mined by the number jn^ of spins that are opposite 
in orientation between basis states \xi) and \xk). For 
example, if |x) = (jlll) -b |lil) + |lll))/-\/3, then 
jik = 2, so that Eq. (2) takes on the form T{t) = 
yJZ -b 2 M[cos 0 i 2 t -b cos 0 i 3 t -b cos023t]/3, where Oik hap¬ 
pens to be 2{B^ — B^). After a semiclassical evaluation of 
the Overhauser field noise, and using the result of T 2 *(l) 
in Eq. (A2), we find M[cos0ikt] = ^ 

g- 8 [t/T,*(i)]=^ so that B{t) « exp{- 8 / 3 [t/r 2 *(l)] 2 } in 
the short time limit. Thus in this example, T 2 *( 3 ) = 


'72'(1)vW 8: with T(t) = exp{-[t/T 2 *(n)] 2 }. 

Examples of multi-spin decoherence .—With our under¬ 
standing of single-spin decoherence, and with a measure 
(fidelity) of the collective decoherence for a multi-spin 
state jx), we are now in position to clarify the scaling 
of the inhomogeneous broadening time Tf{n) in various 
subspaces of the n-spin system. Below we describe the 
results from several representative classes of jx). 

Case A: single product state .—The simplest multi-spin 
state is a single product state (m = 1). The random 
Overhauser field acting on a product state creates a ran¬ 
dom but global phase (relative to when the nuclear reser¬ 
voir is absent). This global phase does not lead to any 
decoherence, as there is no coherence (phase) information 
stored in a product state to begin with. 

Case B: two product states, with m = 2 and k > 1.— 
The simplest multi-spin state that can undergo pure de¬ 
phasing consists of two product states. Here we choose a 
particular class of |x)b = di\b) -\- d 2 \k), with one product 
state being fully polarized \b) = jl)®", while the other 
being from the fc-th subspace with k spins prepared in 
11). The fidelity of such a state is given by 

dfU). I ’ 


T{t) m exp < —A\did 2 fk 


so that 

Ti{n) 1 
Tf{l) 2\did2\Vk' 


( 4 ) 


In this case, dephasing time is inversely proportional to 
the square root of the number of spins prepared as jl) in 
jfc). A special example here is the GHZ state, |x)ghz = 
(|l) 0 ra _|_ for which the two product states 

have completely opposite spins. The decoherence rate is 
simply Tf{n)/Tf(l) = where n is the number of 

spin qubits involved. Indeed, the worst case of scenario 
for a two-product-state |x) is when the pair of product 
states have completely opposite spins, Tf{n)/Tf{l) = 
l/yfd. 

Case C: n > m > 2,k = 1—We now consider an jx) 
that is a general superposition of m product states from 
the first manifold with one spin in jl). In other words, 
|x)c = dl|lll---I) + d 2 |ll---n) + ••• + dnllll---), 
where MjP = 1 - This is a state that is slightly 

more general than the W-state, with a random weight 
and phase for each basis state. The fidelity of |x(t))c' is 


T{t) ft! exp 


® 'y Mil *^121 

ii <h 



( 5 ) 


which implies (by the Cauchy—Schwarz inequality) 


1 I n Tf{n) 
2y n-1- Tf{l) 
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< 00 . 


( 6 ) 


Here the upper bound (00 means no decoherence) is ap¬ 
proached when a particular \dj., \ = 1 while all other 
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so that we go back to a single product 
state. The lower bound corresponds to the equally- 
populated superposed states with \dj\‘^ = 1/n, i.e., an 
almost normal IT-state (a standard W state would have 
all dj having the same phase, too). When n —>■ oo, 
r 2 *(n)/T 2 *(l) > 1/2. The whole system acts like a gi¬ 
ant spin—1/2 system that is spread out over n physical 
spins. Notice that the lower bound of decoherence time 
is scale-free. The scaling of decoherence for large n is 
insensitive to either the population distribution on each 
basis state or the total number of physical spins. 

Case D: m = C^,k> 2.—We now extend |a;) to a more 
generalized W-state that is uniformly distributed over all 
the product bases in the /c-th Zeeman manifold, with k > 
2. Consider a particular example, \x)d = , 

where \dr\^ = l/C^” and each \xr) has k spins in |1) and 
n — k spins in |1). The overall decoherence is determined 
by the phase differences between every pair of states from 
the basis states. Since we limit our 

discussion below to k < nj^ without loss of generality. 
The phase difference between a particular pair of 

product states Ix^) and 1x^2) can involve Overhauser 
fields in 2j QDs, where / < k. In the extreme case of 
2j = n, the pair of states has completely opposite spins. 
For each possible j, there are pairs of states 

with the same phase difference as well as the same value 
of ensemble average = exp{— 2 j[ 2 t/T 2 *(l)]^}. 

Therefore the fidelity for this generalized W-state is 


Ht) = + E 

” y ri<r2 
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exp ■ 


t 


Thus 


Akin — k) 

T^{n) 1 


T*{1) 2\^k(n-k) 


(7) 


( 8 ) 


In this case, we find that (a) when n —^ oo, T^iji) « 
l/(2-\/fc), which is scale-free with respect to the number 
of spins n as well as the number of product states m in 
\x)d (it is a similar feature as in Case C, where fc = I); 
(b) overall decoherence is completely suppressed when 
fc = 0 or fc = n, i.e. T^iji = k) = oo. These two Zee- 
man manifolds contain one state each, so that Case D 
is reduced to Case A; (c) the strongest decoherence oc¬ 
curs when k ~ n/2, where T^in) « T|(I)/-yn; (d) the 
generalized VF-state |x)d here is a reliable lower bound 
for the decoherence scaling rate of a more general state 
\x)'jj = 'Y^^=i^r\xr) in the k-th manifold. In Fig. 2, 
the dashed lines represent the analytical result given in 
Eq. (8) with 7 < n < 20 and 2 < fc < 4, and the data 
for error bars are obtained by the maximal errors numer¬ 
ically evaluated on 100 IxYd states randomly generated 
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FIG. 2. The result of T 2 {n)/T 2 (1) by randomly generated 
\x)'o states in the fc-th manifold comparing to the analysis 
expression by the generalized TF-state in Case D. As com¬ 
pared to the equally-populated |x )d state, a |x)/) has random 
populations in each product basis state. 



in the fc-th manifold, with completely random coefficient 
for each basis state. Figure 2 clearly shows that devia¬ 
tions from the result of |x)D-state quickly decrease with 
increasing n and k. Thus the equal-weight |x)d state is 
a very good representative of both Cases C and D. 

Case E: m = n .—Next we further generalize |x) to 
be a superposition over product states picked from more 
than one Zeeman manifold. Out of the infinite number of 
possible combinations, we pick one class of such states, 
with one product state from each Zeeman manifold, so 
that m = n and |x)£; = J2r=i^r\^r)j with \xr) picked 
from the r-th manifold. To obtain analytical results, 
we first assume equal weight for all the states involved: 
\drY = ^/n. To further limit the choice of states, we as¬ 
sume there is j spin polarization difference between \xr) 
and |xr+j), r = l,2,-- - ,n — j. For example, for a 3-spin 
system, |x) can be chosen as (|111) -I- |ill) -|- |ni))/-\/3. 
For an arbitrary n, the fidelity is 


Tit) = 




exp < -- 


- j) exp < - j 
i=i I 

2(n2 l) 


2t 


lTUl)\ 


3n L7^2(1)J 


(9) 


Therefore 


Ti{n) I 3n fY 

T^{1) ~ Y 2(n2 - I) V ^ 


( 10 ) 


As in Case D, we generalize \x)e to |x)^ by random¬ 
izing the weight |dr P’s, I < r < m, and the selection of 
the product basis states within each manifold. In Fig. 3 
we plot our numerical results as compared with the ana¬ 
lytical expression from Eq. (10). While the error bars in 
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FIG. 3. The results of TJ (’i-)/T '2 (1) by randomly generated 
\x)'e states over all the Zeeman manifolds, as compared with 
Eq. (10) in Case E. In frame (a), we use 100 states with 
random product bases and random populations; in frame (b), 
we use 100 states with the fixed bases of \x)e but random 
populations in each product basis state. 
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Fig. 3 for random are larger than those in Fig. 2, 

the analytical result is still a good indicator of the aver¬ 
age T^in). The size of the error bars also decreases with 
increasing n. Thus the n-spin dephasing time scales as 
l/v^ for a large n. 

Case F: m = 2" — 1.—Suppose |x) is a combination of 
Cases D and E: |x)f = \/C'^/2”|Wfe), where \ Wk) 

is a normalized bF-state in the fc-th manifold. In fact, 
|x)f is just the fully superposed state [(|1) + |l))/-\/2]®". 
For the overall decoherence. Cl™ pairs of phase differences 
have to be taken into account. There are 2"C'^ elements 
involved with j spins, 1 < j < n, in the set of *i 

i 2 . The fidelity is 


F{t) Ri exp 



Consequently, we have 

T^jn) 1 

r*(i) ^ 


( 11 ) 


( 12 ) 


Based on this equation and the numerical simulation over 
|x)^ with randomized coefficients as seen in Fig. 4, the 
dephasing time for |x) residing in the whole Hilbert space 
adhere to the sublinear power-law the same as in 

Case E. 

Conclusions and Discussions .—We have explored the 
scaling behavior of the decoherence time of n decou¬ 
pled electron spin qubits by investigating the fidelity of 
6 classes of representative superposed states |x). Each 
electron spin is individually coupled with its own nuclear 
spin bath through hyperfine interaction, and we do not 
consider electron-electron interactions in this study. 

Our results are summarized in Table I, where k is the 
number of spins in |1) in a product state that makes 


lx) 

n{n)/n{l) or r 2 (n)/r 2 (l) 

Stable: A 

no decoherence 

Two product states: B 

fc-5 

fc-th subspace: C and D 

fc-5 

Crossing subspaces: E and F 

I 

n 2 


TABLE I. Various scaling behaviors of decoherence time of 
the n decoupled-electron-spin system under the influence of 
hyperfine coupling with local nuclear baths. 


up of |x). Typically, both inhomogeneous broadening 
dephasing rate l/TJ"(n) and pure dephasing rate 1 /T 2 (ri) 
are sublinear power-law functions of spin number n. If 
|x) is constrained in a single subspace with a fixed fc, 
T^ln) and T 2 {n) become scale-free with respect to n and 
m (the number of basis states involved). 

The scaling behaviors revealed in our case studies can 
be qualitatively understood based on counting the num¬ 
ber of different spin orientations in any pair of product 
states. Considering any m product states making up a 
|x), a large fraction of pairs have 0{n) electron spins ori¬ 
ented in the opposite direction. If we average over all 
possible states assuming jdrp ~ I/w, then the state fi¬ 
delity given in Eq. (2) could be estimated as 


F{t) 


1 m{m 1) 


exp 


m — 1 
m 


0{n) 



The decoherence rates are insensitive to m because 
of normalization and our equal-population assumption. 
Furthermore, in the fc-th manifold, the scaling law is 
1 /Vk because an arbitrary pair of states is different in 
0(fc) spins. 

Our study here could be straightforwardly extended to 
other decoherence mechanisms. If the single-spin deco¬ 
herence function is given by W{t) = exp{—[t/T 2 (l )]*"}, 
the index of every power-law (—1/2) in Table I should be 
revised to —\/v. For example, spin relaxation induced 
by electron-phonon interaction produces a linear expo¬ 
nential decay characterized by Ti, with I/T 2 = l/(2Ti). 
In this case the scaling power-laws for the n-spin system 
will be modified to be proportional to k~^ or n~^ based 
on the selection of |x). For decoherence due to Gaus¬ 
sian noise under dynamical decoupling [27], the decay 
functions have z/ = 4 for spin echo (SE) and r/ = 6 for 
two-pulse Carr-Purcell-Meiboom-Gill sequence, so that 
the scaling factors for decoherence times of the n-spin 
system become and respectively. 

Our results are important to the scale-up considera¬ 
tions for spin-based quantum computers or more general 
qubits that are under the influence of local reservoirs. 
The sublinear scaling shows that a large superposed state 
does not lose its fidelity overly quickly as conventional 
wisdom may dictate. The scale-free states also help us 
identify what Hilbert subspaces are more favorable in co¬ 
herence preservation. 
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Appendix A: Single-spin Decoherence 


For a single electron spin coupled to the surrounding 
nuclear spins in a finite magnetic field, the nuclear reser¬ 
voir causes pure dephasing via the effective Hamiltonian 
[16, 17] 


= 2S^{Ha + U) 


N 


= 2S- ^ 


An 




E 

ol^ol' 


AoiAfy 


-^^4' -(Al) 


where N is the number of nuclear spins, H is the elec¬ 
tron Zeeman splitting, and is the hyperfine coupling 
strength. The sums over a and o! here are over all the 
nuclear spins in the single quantum dot (QD). The de¬ 
phasing dynamics has two contributions: Ha is the lon¬ 
gitudinal Overhauser field, while V is the second-order 
contribution from the transverse Overhauser Held. In a 
finite field, normally the former dominates, generating a 
random effective magnetic field of AB ~ 1 to 5 mT [12] 
on a quantum-dot-confined electron spin in GaAs. This 
random field leads to a stochastic phase and accounts for 
the inhomogeneous broadening effect characterized by a 
free induction decay at the time scale of T 2 *(n), with n 
being the number of single-electron QDs. For a single 
dot n = 1, the inhomogeneous broadening decoherence 
function is: 


M 




= W^(t) = exp 


L 441)4 


■ (A2) 


Appendix B: Notations on the Overhauser fields 


A convenient way to understand the effect of hyper¬ 
fine interaction on the n-decoupled-qubit system [see 
Eq. (1)] is to introduce the semiclassical Overhauser field: 

~ where p = z,+,— refers to the 

longitudinal and transverse directions, Ij takes the value 
of 1 or 1, and Bj = a Overhauser field 

in the jth QD. In a finite field and up to second order, 
the hyperfine Hamiltonian could be diagonalized on the 
product state basis into 


Hm 


-diag 






(HI) 


where 


B 


Z 

— 1 ‘ ‘ ‘^1 


= B 


z 

l-nln — 1" 'll 


1=1 


(B2) 


Here pj = + {—) ii Ij = 1 (I). The two terms in 
Eq. (B2) are responsible for the inhomogeneous broad¬ 
ening and narrowed-state FID, respectively. Accurate to 
the first order in Aja, B~B^ fv B^B~ since [B^ ,B~] = 
2 is second order in the hyperfine coupling 

strength and is small. For simplicity we take pj = + 
in the following derivation. Generally, the correction 
term for the jth dot in Eq. (B2) Ri ljB^B~/(2fl). For 
example, a completely polarized state \b) = |I)®" ex¬ 
periences a longitudinal Overhauser field B^ = B§ + 

EUb+b-/{ 2n). 

Now that the hyperfine Hamiltonian takes on a diago¬ 
nal form, it can only lead to dephasing between different 
product states due to B, similar to the single-spin case we 
discussed above. The dephasing of a product state [xr) 
relative to \xr') is due to the difference in the random 
Overhauser field B for these states. 


Here M[-] is an ensemble average over the longitudinal 
Overhauser field in the QD, and T^il) oc '/N/A with 
A = Aa- In a single gated QD in GaAs, T^ll) is in 
the order of 10 ns. 

If the effect of Ba is suppressed, such as through nu¬ 
clear spin pumping and polarization [13], V, which is sec¬ 
ond order in the transverse Overhauser field, leads to the 
so-called narrowed-state free induction decay (FID), by 
which the off-diagonal elements of the spin density ma¬ 
trix decay at the time scale of In the main text 

we simplify the notation for to T 2 {n), in the same 

way as T^in). The narrowed-state decoherence function 
for a single dot is: 


W^{t) 


exp 


e ■ 

4(1)4 


(A3) 


Appendix C: Statistical independence of 
inhomogeneons broadening and narrowed-state free 
induction decay 

To analyze the relationship between inhomogeneous 
broadening and narrowed-state free induction decay in 
an n-decoupled-qubit system, we consider an arbitrary 
pure state in a subspace spanned by m spin product 
states Jx) = J2T=i^r\xr), where [xQ = 

Here Ij refers to the electron spin orientation along the 
z-direction in the jth QD for state jxr), and takes the 
value of 1 or I = — 1 for notational simplicity. This selec¬ 
tion is general enough to cover all the cases discussed in 
the main text. Helped by the Overhauser fields defined 
above, and under the diagonalized hyperfine interaction 
Hamiltonian in Eq. (Bl), an initial state jx) evolves into 


m 

where 72(1) ^ Nfl/A^ [17], and is in the order of ps in 
a gated GaAs QD. 


\x{t)) = '^dre (Cl) 
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where Collective decoherence emerges 

due to the non-stationary random phase differences 
among the m product states |a;r)’s. The fidelity between 
|x(0)) and |x(t)) can be expressed as 


The ensemble average M[cos9kri\ = could be 

estimated using the decoherence times of a single qubit 
system T 2 *(l) (inhomogeneous broadening time scale) 
and T2(1) (the narrowed-state FID time scale) [17], 


T{t) = \/M[(x(0)|a;(t))(x(f)|a:(0))] 


(C2) 


\ 


M 






\r—l 


\r—l 


\ 


M 


E ^E cos Okrt 


r—1 


k<r 


(C3) 


where the phase differences 9kr = — B^^. Accord¬ 

ing to Eq. (B2), each 9kr could be decomposed into two 
terms, 9]^^ and 0^®, that are responsible for the inhomoge¬ 
neous broadening and narrow-state free induction decay, 
respectively: 


n _ nib , nns 

^kr — ^kr "T ^kr 5 

nib _ TDZ 

^kr ~ 

1 




n ‘'n — 1 *'1 


= E(7 - 'PSJ. 




^kr 


2D 




i=i 


M[cos 9 krt] = M 


e<rte<rt 


= M 


e 


M 


e 


i=i 


t2 

+ 


T2*(1)2 T2(1)2 


.(C4) 


This result is obtained by the quantum noise theory [28] , 
which is valid at least in the short time limit. Physi¬ 
cally it is based on the assumption that longitudinal and 
transverse Overhauser fields are independent from each 
other, so that the averages above can be factored. The 
two decoherence mechanisms are thus mutually indepen¬ 
dent. Using the short notations Bkr = ~ 

= [t/T^*{l)]^ and B"®(t) = [t/T 2 {l)]^ Eq. (C3) 
can be rewritten as 


Bit) 


Y, \dr\* + 2^ 14Pld,|2 [1 - BkrD-b^t) - Bfe,B"«(f)] 

\ r-1 


k<.r 


1- hYldkl^ldrl^Bki] [D^^it) + D"«(t)] = - B[B*pt) + D-®(t)] 

\ \ k<r / 


exp 




exp 


' B 


(C5) 


where B = 2j2k<r\dk\^\dr\'^Bkr- In short, Eqs. (C4) 
and (C5) prove that inhomogeneous broadening and 
narrowed-state FID are independent decoherence chan¬ 
nels, and have the same scaling behavior. The over¬ 
all decoherence function is just a simple product of the 
decay functions for inhomogeneous broadening FID and 
narrowed-state FID. We can thus focus on just inhomoge¬ 
neous broadening in our discussion of decoherence scaling 
for n spin qubits and in main text, we omit the super¬ 
script of the phase difference 9 for notation simplicity. 


Appendix D: Numerical evaluation of n-spin 
decoherence 

Equation (C5) gives a general description of decoher¬ 
ence function within the Overhauser field approach. No¬ 
tice that the fidelity Bit) of a pure state |a:) is solely 


dependent on the function B, which is only a function of 
populations in product basis |a;,.)’s, but not a function of 
the phases of these amplitudes. For example, when jcc) 
is a single product state |xr), i.e., Case A in the main 
text, B = Q. Now the random phase is global, and does 
not lead to decoherence. If more than one coefficient is 
non-vanishing, so that B ^ 0, there will be finite deco¬ 
herence. 

One can then use the expression of B and Eq. (C5) to 
numerically obtain the scaling behavior of T 2 *(n)/T 2 *(l) 
or T' 2 (n)/r 2 (l) for an arbitrary initial state. Although 
there is an infinite number of possible superposed states 
even for a small n, it turns out that the averages of the 
numerical results for different classes of states agree quite 
well with the analytical expressions obtained in the main 
text. For example, in Fig. (4), we use 100 randomly 
generated states with the bases of \x) in Case F in the 
main text, i.e., all the product basis states in the whole 
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Hilbert space, but random populations in each product 
basis state. The error bars in Fig. 4 for random states 
rapidly vanishes with increasing n. Making the analytical 
expression a really good predictor of decoherence for an 
arbitrary state. 


FIG. 4. The results of T 2 {n)/T 2 (1) by randomly generated 
states over the whole Hilbert space of n-spin system, as com¬ 
pared with Eq. (12) in Case Fin the main text. For each n, 
we use 100 randomly populated states. 
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